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Abstract 

We study the transverse momentum dependent factorization for single spin asymmetries in Drell- 
Yan and semi- inclusive deep inelastic scattering processes at one- loop order. The next-to- leading 
order hard factors are calculated in the Ji-Ma-Yuan factorization scheme. We further derive the 
QCD resummation formalisms for these observables following the Collins-Soper-Sterman method. 
The results are expressed in terms of the collinear correlation functions from initial and/or final 
state hadrons coupled with the Sudakov form factor containing all order soft-gluon resummation 
effects. The scheme independent coefficients are calculated up to one-loop order. 
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Single transverse spin asymmetry in hadronic reactions have attracted great attentions 
from both experiment and theory sides in recent years. It promises strong connection to the 
three-dimensional partonic tomography of the nucleon, and provides unique opportunities 
to study QCD dynamics, such as the factorization and universality of parton distributions. 
In particular, the short distance partonic interactions resulting into different Wilson lines 
in the relevant transverse momentum dependent (TMD) parton distributions between the 
Drell-Yan lepton pair production in pp collisions and the semi-inclusive hadron production in 
deep inelastic scattering (SIDIS) predicts opposite signs for the Sivers-type single transverse 
spin asymmetries in these two processes [H, 0]. This nontrivial universality prediction has 
stimulated strong interests around the world to measure, especially, the SSA in the Drell-Yan 
(DY) process, since that in the SIDIS has been observed in various experiments 0]. 

However, the numerical predictions for the SSAs in these processes are all based on 
a leading order naive TMD factorization Q. Although a general TMD factorization has 
been argued for these processes \^$\, we need to know the next-to-leading-order (NLO) 
perburbative QCD corrections to have more reliable predictions, which have not yet been 
calculated. Another important issue is the energy dependence. Current DIS experiments 
cover the Q 2 range about 2-5GeV 2 , whereas that in the planned DY process will be measured 
at relative larger Q 2 about 20-25GeV 2 . Here, Q represents the large momentum scale, i.e., 
the virtuality of the photon in these processes. To accurately describe the Q 2 evolution of 
the transverse momentum dependent observables, the QCD resummation effects have to be 
taken into account (6|. 

In this paper, we will build a theoretical framework to address these important ques- 
tions. We carry out, at the first time, the complete NLO perturbative correction to the 
single transverse spin dependent cross sections in DY and SIDIS processes in the TMD fac- 
torization. One of the important implications of the explicit one-loop calculations is to help 
to construct the correct resummation formalism for the SSA observables. Earlier attempts 
to formulate these effects in SSAs have been made in various forms 111! . A resummation 



formula close to ours was used in Ref. [10j . and a significant suppression effect were found 
when Q 2 is very large. In the following, based on the one-loop calculation results, we will 
derive the complete soft gluon resummation formalism following the Collins-Soper-Sterman 
(CSS) method Q, which can be easily implemented in the phenomenological studies. 

We take the DY process as an example to demonstrate our procedure and present the 
main results, 

a(p a , s±) + b(p b ) 7*0?) + x t + r + X, (1) 

where Pa and Pb represent the momenta of hadrons A and B, and S± for the transverse 
polarization vector of A, respectively. The single transverse spin dependent differential cross 
section can be expressed as 

^L = ^ sl wUQ^), (2) 

where q± and y are transverse momentum and rapidity of the lepton pair, o"o = A7ral m /3N c sQ 2 
with s = (Pa + Pb) 2 , and e al3 is defined as e a ^^ v Pa^Pbv / Pa • Pb- At low transverse mo- 
mentum (q± <C Q) the structure function Wjjt can be formulated in terms of the TMD 
factorization where the quark Sivers function is involved [l|, [2[ , whereas at large transverse 
momentum (q± 3> Aqcd ) it can be calculated in the collinear factorization approach in 
terms of the twist-three quark-gluon-quark correlation functions 12|~14]. It has been shown 
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that the TMD and collinear twist-three approaches give the consistent results in the inter- 



mediate transverse momentum region: Aqcd "C q± <C Q [14], [15] • This consistency allows 
us to separate Wut into two terms [6|, 

/i 2 i 
—-e^WS T (Q; b) + Y&iQ; q ± ) , 

where the first term dominates at q± <C Q region, and the second term at q± ~ Q. The 
latter is obtained by subtracting the low q^ expansion from the full perturbative calculation. 
According to the TMD factorization, we have 0, EH, 

WdAQ;b) = f£ a \ Zl ,bXi)q(z 2 ,bX2) 

xHur^iS^p))- 1 , (3) 

where z± t 2 = Q/y/se ±y and the sums over flavor weighting with the charge squared of the 
quarks are implicit, and q are the TMD quark Sivers function of A and antiquark 
distribution of B, H and S are hard and soft factors, respectively. In this paper, we follow 
the Ji-Ma-Yuan factorization scheme 0], where two off-light-front Wilson lines (along off- 
light-front vectors v i and v 2 ) are introduced to regulate the light-cone singularities for the 
TMD quark distributions. We further define (f = {y\ ■ Pa) 2 jv\ and ( 2 = {v 2 ■ Pb) 2 /v 2 , and 
the rapidity cut-off parameter p: zf(f = z%C 2 = pQ 2 . The Sivers function in the impact 
parameter 6^-space is defined as (x, b i ) = j d 2 k±e~ t ^ ± ^ ± k^_f^ DY \x, k±) /Mp, where 
fyr follows the definition of Ref. [16J. Our results can be translated to other factorization 
schemes where different regularizations for the light-cone singularities are used j9|. 

We investigate the above factorization formula Eq. (3) in the perturbative region of 
> Aqcd- The explicit calculations at one- loop order will verify the TMD factorization, 
from which we obtain the NLO correction to the hard factor. In the calculation of the 
spin-average (or double spin asymmetry) TMD factorization, it is convenient to choose a 
quark (or gluon) target, where every factor in the factorization formula can be calculated 
perturbatively (7|. However, the SSA vanishes with on-shell quark. To get nonzero effect, 
we have to go beyond the simple quark target picture. 

Our calculations are based on the collinear correlation functions from the incoming 
hadrons. In this framework, the SSA is naturally a twist-three effect, and involves the twist- 



three quark-gluon-quark correlation function from the polarized nucleon [12|, |13J . There have 
been great theoretical developments in this framework in recent years, see, e.g., Refs. (ljini- 
20|. We will utilize these techniques to compute the structure function W{j T (Q,b). First, 



let us write down a general form, 

W{; T (Qi fy = —± / 2 , T F (x 1 ,x 2 )q(x > ) 

V 2 J J x 1 x 2 x' 

xH(x 1 ,x 2 ,x';Q,b) , (4) 

where Tp is the Qiu-Sterman matrix element for quark q 13, 1^| and q(x') the integrated 
anti-quark distribution. The quark Sivers function is related to the Qiu-Sterman matrix 
element: Tp(x,x) = J d 2 k±\kWf^ DY \x,k±)/Mp [H], 19]. Other twist-three quark-gluon- 
quark correlation functions contributing to the SSA can be included as well. For simplicity, 
we focus on Tp contributions in this paper. 

The above formula is similar to that of the q± -weighted asymmetry in the same process 



calculated in Ref. [191 ] . There are two major differences: (1) instead of weighting with g", 
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FIG. 1. Generic Feynman diagram contribution to the impact parameter space structure function 
W UT (Q,b). 



here we weight with e~ iq± ' bj -; (2) we will restrict our calculations of the differential cross 
section in the TMD domain, i.e., q± C Q (V&± *C Q)- We calculate Eq.(4) in a covariant 
gauge, where a generic diagram is illustrated in Fig. 1. Because of the twist-three effect, 
an additional gluon attachment has to be taken into account for hard partonic scattering 
amplitude. A collinear expansion will be performed to calculate the hard part H(Q,b). 
In particular, the amplitude is expanded in terms of kg ± = k% ± — k"j_ where ki± <C l/b±. 
Combining this expansion with matrix element from the polarized nucleon, it will lead to 
the spin dependent cross section expressed in terms of T F (x 1 ,x 2 ). 

The leading order diagrams are shown in Fig. 2(a) and (b). From the kinematics, we 
find that qj_ is related to the transverse momenta of the two quark lines as: q± = h%± for 
Fig. 2(a) and q± = ki± for Fig. 2(b). Therefore, the contributions from these two diagrams 
will be 



Fig.2 
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— e 



-iki 



(5) 



where the last equation comes from the collinear expansion of the exponential factor. The 
initial state interactions represented by the propagator leads to a pole contribution. After 
taking the pole, we will obtain the leading order contribution, 



W^\Q,b) 



T F (z 1 , z 1 )q(z 2 ) 



This also normalizes the leading order hard factor as = 1 in Eq. (3), because 



(6) 



la(DY) ( 
IT 



Zi,b±) = Tp(zi, z%)(— ibj_/2) at this order at small bj 



Order a s corrections come from real and virtual gluon radiation contributions. Similar 
to that calculated in Ref. [19], the virtual diagrams contain soft and collinear divergences. 
The soft divergence will be cancelled out by the real gluon radiation contributions, and 
the collinear divergences will be absorbed into the parton distributions in both TMD and 
collinear factorizations in Eqs. (3) and (4). The real diagrams are similar to those calculated 



in Refs. [14J, [19(. In order to obtain their contributions to Wfi T (b±), we need to first take 
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FIG. 2. Leading order Born diagram calculation of Wut(Q,°). 



the low transverse momentum limit, and then Fourier transform into the impact parameter 
6j_-space. After summing up both real and virtual diagrams contributions, we indeed find 
that the soft divergences cancel out between each other, and the total contribution only 
contains the collinear divergences, 
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23 , 



(7) 



where 6 — z il z [-> Pgg * s the quark _splitting function, and Vqg^ qg the splitting function for 
the Qiu-Sterman matrix element (Til . 2p| . After subtracting the collinear divergences from 
the splitting of Tp and q, we demonstrate the factorization form in Eq. (4) up to one-loop 
order. 

The above result can also be casted into the TMD factorization formula Eq. (3), where 
we have to subtract the TMD quark Sivers function, antiquark distribution and the soft 
factor. The latter two have been calculated before |7|. The quark Sivers function can be 
calculated similarly, and we find that, 



fi T (z 1 , 6j 



2tt 



-ibl 



qg^gg 
2 



T F (^) + 5(l-6)Cf 



'2 ln 6V' 



-27B 



zfCfb 2 



3 27S-1 



7T 



+ 



2N r 



(1-6 



(8) 



After subtracting these factors out, we find that the hard factor Hut is free of infrared 
divergence, and it is the same as that for the spin average one calculated in Ref. Q, 
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which is very interesting and suggests that the hard factors are spin- independent. Following 
the same procedure, we find that the hard factor for the single transverse spin depen- 
dent cross section in SIDIS. Again, it is the same as the spin-average case, ^utIsjdis = 
-f^jylsiDis 01 • Clearly, the hard factors for DY and SIDIS differ by a factor of it 2 . The hard 
factors in other TMD factorization scheme can be calculated similarly. 

The factorization formula Eq. (3) contains large logarithms in terms of In 2 (Q 2 b 2 ) [ll[ 
which is also shown in the total contribution result of Eq. (7). To resum these large loga- 
rithms, we follow the CSS method 0], 



W$ T (Q;b) = e-^^WSAC^b) 
= (-ib°[/2) e-^WEij 

xACj ® fl%fa)Cu ® f 3/3^2) , (9) 

(3) 

where fj/B represents the integrated parton distribution from hadron B, and the twist- 
three function from hadron A (Tp(zi,z 2 ) is the most relevant one). The last step of the 
above equation comes from further applying the collinear factorization formula Eq. (4) at 
lower energy scale C\jb. From the above one- loop calculations, we find that the perturbative 
Sudakov factor Sjjt have the same form as that for the spin-average case, 



S UT {Q\b) 



2 y 



>2 



dfi 2 



In ( ) A UT {C x] g^)) 



+ B UT (C 1 ,C 2 ;g(n))} , (10) 



up to one-loop order, where G\ and C 2 are constants in the order of 1. 

The A, B and C functions can be calculated in perturbation theory: A = J2 n =i A^ n \a s /Tr) r 
From the explicit one-loop calculations, we obtain the following results for these coefficients, 

A™ = C F , £#t = -3/2C7p, AC?® = 6(1 - x) , 



AC2 1 ) = -4r(l-z) + — S(x 
AN r ' 2 v 



^-4 
2 



where C qq follows the spin-average case |6j, and we have chosen the canonical values for 
C\ = 2e~ lE and C 2 — 1 to simplify the above expressions. For SIDIS, A, B remain the 
same, whereas AC T have opposite sign and there is no ir 2 term in AC^ . The above 
coefficients can also be calculated by comparing the fixed order calculations of the differential 



cross section depending on transverse momentum [I4| to the expansion of the resummation 
formula Eq. (9). We have checked that this gives the consistent results. 

Eq.(9) is our main result for the QCD resummation for the single spin asymmetry in the 
DY process. The structure, in particular, the pre-factor is the unique feature for the 
single transverse spin azimuthal angel dependent cross section. This pre-factor comes from 
the explicit one-loop calculation above. It also guarantees the proper behaviors for the spin 
asymmetry at small and large transverse momentum. At very small q±, the asymmetry has 
to vanish, so that the differential cross section and W{j T will be proportional to q°[. On the 
other hand, at large q±, the asymmetry is power suppressed by l/q± from the perturbative 
calculation, which leads to W$ T oc q±/q±- With the correct coefficients extracted from our 
one- loop calculation, these behaviors are satisfied in Eq. (9). We would like to emphasize 
that the TMD and collinear factorizations for Wff T (B) in Eqs. (3) and (4) are crucial to 
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obtain the final resummation formula of Eq. (9). Without these factorization results, we 
can not apply the CSS resummation. 

Another important feature is that Eq. (9) depends on the integrated parton distributions, 
which are universal. The opposite sign for the SSAs in DY and SIDIS is reflected by the 
opposite C coefficients in this formula. Moreover, our resummation formula is scheme- 
independent, although the TMD factorization Eq. (3) depends on the scheme of how to 
regulate the light-cone singularities in the TMD distributions. This can be clearly seen from 
the disappearance of p and Q in Eq. (9) with he above coefficients. 

In summary, we have derived the CSS resummation formalism for the single spin asym- 
metries in DY and SIDIS processes. The relevant coefficients are calculated up to one-loop 
order. These results shall be further studied to understand the energy dependence of the 
SSAs in these processes, and provide more accurate predictions for the DY process which 
is actively pursued by several experiments. Our results shall shed light on all other fc^-odd 
observables, and should be applied to the azimuthal angular dependent observables in DY, 
SIDIS, and e + e~ annihilation processes. We performed our calculations in the framework 
of the collinear correlation functions of hadrons. This allows us to compute the cross sec- 
tions in a consistent and rigorous way. We noticed that recently, a different framework has 
been developed where a gluonic degree of freedom is included in the leading order quark 
target (2l|. It will be interesting to apply this method and compare with our results. 
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DE-AC02-05CH11231. We are grateful to RIKEN, Brookhaven National Laboratory, and 
the U.S. Department of Energy (Contract No. DE-AC02-98CH10886) for supporting this 
work. 
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